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Abstract 

We study low-lying excitations in the ID S = 1 antiferromagnetic valence- 
bond-solid (VBS) model. In a numerical calculation on finite systems the 
lowest excitations are found to form a discrete triplet branch, separated from 
the higher-lying continuum. The dispersion of these triplet excitations can 
be satisfactorily reproduced by assuming approximate wave functions. These 
wave functions are shown to correspond to moving hidden domain walls, i.e. 

to one-soliton excitations. 
PACS numbers: 75.10. Jm 
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I. INTRODUCTION 



It was almost a decade ago that Haldane |J conjectured the existence of a new type 
of ground state for isotropic Heisenberg antiferromagnets (HAF) of integer spin S. The 
Haldane phase was proposed to be characterized by a unique, disordered ground state with 
exponential decay of the correlation functions and a finite energy gap to the excited states. 
Since then one-dimensional quantum spin chains with S = 1 have been studied intensively 
and it is claimed that the gapful behavior is a generic feature of integer-spin models @. 

The first rigorous example of a S — 1 antiferromagnetic model with Haldane phase was 
given by Affleck et al. [|J. They showed that the S = 1 isotropic bilinear-biquadratic model 
defined by the Hamiltonian 
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has a short range valence-bond- solid ground state for (3 = 1/3 (VBS or AKLT model). At 



that point ^(1/3) is a special projection operator, that projects out the quintuplet state of 
the two neighboring spins, and it is positive semi-definite. Therefore, a state Q for which 
hjfl = for any j is necessarily a ground state with ground-state energy E GS = 0. Such 
an Q state could be constructed using nearest-neighbor valence bonds. They were also able 
to prove rigorously that in the infinite chain limit this state is the only ground state, it 
is separated by a finite gap from the excited states and the two-point correlation functions 
decay exponentially. 

According to Haldane's conjecture such a phase should not appear for half-integer values 
of the spin. This was proven rigorously by Affleck and Lieb |4| and independently by Kolb 
H . It was shown for a wide class of models that in the case when the ground state is a spin 
singlet, the energy spectrum as a function of momentum k is symmetric under reflections 
with respect to k = ran/ 2 (m integer), and therefore the ground state should be at least 
doubly degenerate. 

A similar proof fails in the integer S case, allowing for the existence of a unique singlet 
ground state f4j,||. The excitation spectrum is in general symmetric with respect to k = 



ititc only. Of course, higher symmetry can also appear in integer S chains, as e.g. in the 
spontaneously dimerized phase of the general bilinear-biquadratic S — 1 model |J. 

The S = 1 bilinear-biquadratic model is integrable [[!]] at the critical point (3 = 1, that 
separates the dimerized phase from the Haldane phase. Although the spin is integer, at this 
point the excitations can be described in exactly the same way as for the spin- 1/2 HAF. 
More generally it is known since the work by Faddeev and Takhtajan [fjj, that there are 
integrable spin models for arbitrary S in which the elementary excitations are in fact spin- 
1/2 solitons with a dispersion independent of the spin length S. The observable excitations 
are composite particles, since due to topological reasons the solitons can appear in singlet or 
triplet pairs only. As the energy of such a soliton pair can be described by two parameters, 
the excitations form a continuum in momentum space. 

Away from the integrable point, where the symmetry properties of the excitation spec- 
tra are different for integer and half-integer S, the above mentioned picture of composite 
excitations may not hold. In this paper we will study this problem. 

We will restrict ourselves to the S = 1 case only, where the non-integrability appears 
in the most dramatic way in the Haldane phase. We will show that the lowest excitations 
are real spin-1, one-particle excitations, they cannot be decomposed into pairs of spin-1/2 
solitons. These triplet excitations are, however, not usual antiferromagnons but rather some 
sort of hidden spin-1 solitons. 

The solitonic nature of the excitations of the integer spin models was predicted by Hal- 
dane already. This assumption became less astonishing after the discovery of the hidden 
order in the Haldane phase. Recently in a very inspiring work, den Nijs and Rommelse || 
have introduced a nonlocal string operator erg defined by 
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They argued that although in the Haldane phase the ground state is disordered in the 
conventional sense, it has a hidden long-range order that could be characterized by the string 
order parameter 



Osteins^) = . lip h(v%)h , a = x,y } z } (3) 

where h(')h denotes the expectation value in the ground state of the Hamiltonian H. 

This prediction was later verified numerically by several authors ||. The appearance of 
the hidden long-range order was further discussed by Kennedy and Tasaki [TtJ. They showed, 
using a nonlocal unitary transformation, that 0^ ring > corresponds to the spontaneous 
breaking of a hidden Z2 X Z2 symmetry of the model. Similarly, the fact that the four lowest 
states of an open chain are exponentially close to each other is also a consequence of this 
broken symmetry. It is generally expected, that the breaking of a discrete symmetry in the 
ground state leads to an excitation gap, since Goldstone bosons do not appear. Excitations 
of the model can then be thought of as some sort of (hidden) domain walls, separating regions 
with different ground states. This picture was made more explicit by Elstner and Mikeska 



11 1, who used spin- zero defects [12|] to disorder the antiferromagnetic state. The spin- zero 



defects are in fact solitons. One of the main goals of this paper is to further examine this 
problem. 

We will use numerical and analytical methods to study the low-lying excitations in the 
/? = 1/3 case. Beside the fact that the ground state of the VBS model can be constructed 
analytically, there is another good reason to focus on this model. In a recent study of the 
general bilinear-biquadratic model of Eq. ([[]), we observed |§ that the convergence of various 
finite-size estimates to their thermodynamic limit is extremely fast in the close vicinity of 
(3 = 1/3. This is certainly not true for general (3. Moving away from (3 = 1/3 finite-size 
corrections become stronger, and one must consider longer and longer chains in order to 
see the real asymptotic behaviour. The rapid convergence at the VBS point may not be 
very surprising, if we remember, that at this point in the ground state first-neighbor valence 
bonds are only present and the ground-state energy density becomes independent of the 
chain length. Although the excited states do show some dependence on N, this is found 
to be exponentially small for the most relevant levels. Therefore extrapolation from finite- 
size calculations allows to draw quite reliable conclusions on the spectrum and it can be 
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compared directly to the analytical (variational) results. 

After a detailed numerical analysis of the spectrum, we will study trial wave functions 
for the elementary excitations and illustrate their solitonic nature. Since the Haldane phase 
at (3 = is believed to be in the same universality class as the VBS model at (3 = 1/3, our 
finding should be qualitatively correct for the usual HAF. 

The layout of the paper is as follows: In Sec. II we present our numerical results. Section 
III contains the analysis of two seemingly different trial wave functions and their equivalence 
is shown. The elementary excitations are argued to be triplet bonds in the VBS structure. 
In Sec. IV, we recall the nonlocal unitary transformation of Kennedy and Tasaki. The trial 
wave functions for the excitations are studied further, using this transformation, in Sec. 
V. Thus their domain wall nature becomes explicit. Finally, in Sec. VI we summarize our 
results. 

II. NUMERICAL RESULTS 

Our first aim is to study the excitation spectrum of the VBS model numerically, using a 
periodic boundary condition. The symmetry properties of the Hamiltonian in Eq. ([!]) allow 
us to classify the eigenstates according to their total spin St, its component along the z- axis 
S T , and the momentum k = 2irl/N (I integer) of the states, where N is the length of the 
chain. We computed several low-lying eigenvalues of the Hamiltonian for each possible value 
of k, using a Lanczos algorithm, and also determined the total spin of the states. Chains 
with an even number of sites up to N = 16 were considered. 

Figure [1] shows the six lowest eigenvalues for all k in our longest chain with N = 16. For 
some of the energies the total spin quantum number St is also given. It is seen that in full 
agreement with all previous results 0, the lowest excited state (denoted by A) is an St = 1 
state with momentum k = it. Moreover, also in the whole range \k\ > 7r/2 the lowest-energy 
excited states have the same total spin St — 1- In the thermodynamic limit, these states, 
as a function of k, seem to form a continuous branch of excitations. In fact, according to 
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the general theorem we proved in Ref . , an excitation with St > 1 cannot remain isolated 
in the N — > oo limit. 

On the other hand, near k = tt, the energies of the next higher lying excitations are 
situated at a distance from the triplet branch, that is much greater than their average 
energy difference from each other. This behaviour indicates that in the infinite chain limit, 
at least near k = tt, the lowest triplet excitations do not belong to a continuum. The 
existence of such a discrete branch below the higher lying continuum seems to be another 
characteristic feature of the non-integrable integer-spin models. Unlike the integrable spin 
models and general spin models of half-integer S, the lowest-lying excitations in the Haldane 
phase are real one-particle spin-1 excitations. They cannot be decomposed into S = 1/2 
solitons. 

Above the triplet branch, the higher-lying excitations of the VBS model seem to be 
"dense" for all k, supposedly forming a continuum. They probably cannot be described by 
a single parameter. For iV = 16, the triplet branch merges into this continuum somewhere 
below k ~ 7r/2. Near k = there is no sign of a discrete branch, here the lowest-lying 
excitations are thought to belong to the continuum. 

There is another remarkable feature that can be observed in our finite-chain calculation. 
The gap A# to the lowest excited state of the k = subspace is approximately twice the 
singlet-triplet gap A^ at k = tt. The same property was observed by Takahashi |13| for the 
pure Heisenberg chain j3 = 0. Similarly, the gap A c between the ground state and the second 
excited state in the k = tt sector seems to be three times as large as the singlet-triplet gap. 
It is also noteworthy that state B (the lowest k = excited state) is a quintuplet (St = 2) 
state and C (the second lowest k = tt excited state) is a state with St = 3. The physical 
picture behind such a behavior is simple. The excitations near k = can be composed of 
two low-lying excitations near k = tt and similarly, three excitations near k = tt can be 
combined to give another excitation near k = tt. 

Whether the spectrum has this property in the N — > oo limit was tested by extrapolating 
the finite-size calculations to infinitely long chains. In Fig. ^| the gaps Aa, Ab, and Ac are 
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plotted as a function of 1/N. The convergence to the thermodynamic limit is very fast, 
especially for A^. Using standard extrapolation methods, the limiting values of the three 
quantities are: A A = 0.350 124±10~ 6 , A B = 0.71±0.01, and A c = 1.09±0.03, respectively. 
As is seen, the ratios give the anticipated values 1:2:3 within 4% of error. 

This numerical calculation supports rather convincingly the idea that the elementary 
excitations of the model form a discrete triplet branch, which is separated from the multi- 
particle continuum in a wide range around k = tt. Analyzing the lower boundary of this 
multi-particle continuum, it seems very likely that the two-particle states near k = are 
essentially scattering states of two elementary excitations. The energy and momentum of 
such a multi-particle state is then simply the sum of the energies and momenta, respectively, 
of the two particles. The situation is similar for the three-particle states near k = n. In 
the numerical calculation there does not seem to be any sign of bound states below the 
scattering continuum. 

III. TRIAL WAVE FUNCTIONS FOR THE ELEMENTARY EXCITATIONS 

In this Section we will study the elementary excitations of the VBS model analytically. 
We will see that it is possible to reproduce the dispersion relation of the discrete triplet 
branch quite precisely by assuming simple trial wave functions. It will be argued that the 
elementary excitations are (hidden) solitons that destroy the hidden order of the ground 
state. 

First we recall the form of the ground-state wave function. As it was mentioned in 
the Introduction, the ground state of the VBS model can be constructed analytically using 
nearest neighbour valence-bonds. For this an S — 1 operator will be composed of two 
S = 1/2 operators. Taking the tensor product of the two spin-1/2 spaces, a new orthogonal 
basis at site i || is constructed in the form 



where ip a and ipp represent the eigenstates of the two spin-1/2 operators, and the Greek 
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(4) 



indices take the values j and j. There are three independent symmetric combinations 
corresponding to the three eigenstates of the spin-1 operator, |0)», and |— )j with 

Sf = 1,0 and —1, respectively, 

\+)i = if>yV2,] 

|0)i = ^ = ^ T , (5) 

\-)i=i>yv2. 

The fourth orthogonal state that completes the basis is the antisymmetric combination 
corresponding to an S = state at site i. This configuration will be excluded. 

The ground-state wave function of an open chain of length N can be written in terms of 
these states as 

n( ai) M = ^e&^tf^eP** . . . ^ Pm+1 ■ ■ ■ <vftv • ( 6 ) 

Here and in what follows summation is meant over repeated indices. e a ^ is an antisymmetric 
tensor with At both ends of the chain there is a loose spin- 1/2 degree 

of freedom, denoted by a± and (3n- Since both can assume any of the two eigenstates 
independently, the ground state is fourfold degenerate. Three of these states constitute the 
three components of a spin triplet, while the fourth state is a spin singlet. It was shown, 
however, that these four ground states converge to the same infinite volume limit as iV — > oo. 

A unique ground state can be formed even for finite in the case of periodic boundary 
condition by antisymmetrizing the two loose end spins. The ground-state wave function can 
be written as 

n = Q(a 1 ,(3 N )e^ . (7) 

Note that these states are not normalized, ||f2(ai,/?jv)|| 2 = 5^ + ®W anc ^ II^IP = 3 N + 
0(1). 

Q has the interesting property that the configurations appearing in Q look like in the 
conventional S z representation as 
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...0 + 0... 0-0... + 0... 0-0..., (8) 

i.e., each + is followed by a — with an arbitrary number of states in-between and vice 
versa. This is nothing but a dilute spin- 1/2 Neel antiferromagnet, where the 0's represent 
a background and the + and — states denote the two possible degrees of freedom of a spin- 
1/2 particle. As is seen, the hidden Neel order is perfect for the VBS ground state. For 
this case the string order parameter of den Nijs and Rommelse is || 0^ tring (Hy BS ) = |, 
a = x,y,z. Moving away from the VBS point, quantum fluctuations begin to destroy the 
above structure of Q. However, the hidden long-range order, characterized by 0^ ring > 0, 
a — x,y, z, is expected to persist in a wide region, in the whole Haldane phase ||. 

As for the excited states of the VBS model, our knowledge is much less accurate, since 
the eigenf unctions cannot be constructed in a similarly rigorous way. Recently, however, 
two seemingly rather different trial wave functions were proposed to describe elementary 
excitations in the model. Arovas, Auerbach and Haldane ]l4j] proposed the form 

N 

\k) = N- 1 / 2 J2 eikjS j\V), A* = *,+,- (9) 

3=1 

and obtained 

(1 s (k\H YBS \k) 25 + 15cos(fc) 

<k) = (k\k) = 27 ' (10) 

for the dispersion relation of the excitations. As it is seen from Fig. [2], the variational ansatz 
yields an upper bound < ^ = 0.3704, very close to the real excitation gap obtained 
from the finite-size calculation at k = it. 

The dispersion relation Eq. (|T0"D is plotted in Fig. |I]with a dashed line. Comparison with 
the numerical results suggest that the trial wave function \k) gives a reasonable descriptions 
of the elementary excitations not only at k = tt, but in a large region of the Brillouin zone, 
in the range \k\ ^ ir/2. Below that, in the region near k — 0, the two-particle scattering 
continuum dominates the spectrum. 

Looking at the form of the trial wave function in Eq. (^), one is tempted to interpret the 
elementary excitations as magnons. Note, however, that although these "magnons" would 



not destroy a conventional long range order, they do destroy the hidden order. This is 
directly seen if we use e.g. the \i = + component in Eq. (Q). In Sj~Q only those configurations 
appear that look like 

...0-0. ..0 + 0. ..0 + 0. ..0-0. ..0 + 0... , (11) 

i.e., the Nel order of the nonzero components is broken at one point. In the dilute antiferro- 
magnetic picture this is nothing else but a usual antiferromagnetic soliton embedded in the 
background of 0's. Such a soliton, unlike a magnon, destroys the long range order, in this 
case the hidden antiferromagnetic order corresponding to Og tring . 

The solitonic nature of SjQ is much less obvious at first sight, since in this case the 
hidden Neel order of Eq. (|[) seems to remain intact. However, as it will be illustrated in 
Sec. V, now the hidden order in the transverse directions (i.e. Of tring and Os tring ) will be 
destroyed. 

In an alternative approach Knabe ]15| proposed the following construction for the ele- 
mentary excitations. Let us retain the valence-bond structure of Q for every bond except 
between sites j and j + 1, where the two spin- 1/2 degrees of freedom /3j and atj+i are 
symmetrized to form a triplet bond, instead of the antisymmetrized singlet bond. 

To make this construction more explicit, divide the sites of the chain into two sets: 
£ = {1,2, . . ., j} and K = {j + l,j + 2, . . . , N}. We define first the states $f {a u (3 N ), 
with (3j,aj + i =t, I and 1 < j < N, as explicit tensor products of two arbitrary ground 
states, one on L and the other on TZ, respectively, as 

$f ' aj+1 ( ai ,/3 N ) = n c ( ai ,(3 3 )®n n (a J+1 ,(3 N ) . (12) 

Obviously the ground state of the full chain can be obtained by antisymmetrizing with 
respect to <x,- and Pj+i, i.e., by connecting the C and 1Z sides with a singlet valence bond, 

n(a h p N ) = $jV,#v) - $f(ai,/3iv) . (13) 

On the other hand, symmetrization with respect to ctj and (3j + \ defines three new states, 
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^{a u (3 N ) = ^]\a u l3 N ) , (14) 
^( au (3 N ) = ^=[^/( ai ,(3 N ) + ^( ai ,(3 N )} , (15) 
$j( ai ,p N ) = $ l j i (a 1 ,f3 N ) . (16) 

In these states the singlet valence bond between sites j and j + 1 is substituted by a triplet 
bond. This adds an extra spin-1 degree of freedom to the two free spin-1/2 variables in the 
ground states of the open chain. The factor l/y/2 is introduced in <3>° to ensure that the 
states have the same norm, ||<&?(ai, /3jv)|| 2 = j3 N + 0(1), a = +, 0, — , as N — > oo. 

In the case of periodic boundary condition, the two loose spin-1/2 degrees of freedom at 
the chain ends should again be contracted with an e tensor. This defines the functions 

$!* = ^(a 1 ,f3 N )e ()N ' ai , a = +,0,-, (17) 

for 1 < j < N — 1. The state $^r, in which the spin-1 bond connects the last and first 
sites, N and 1, is defined analogously. Unlike the case of open boundary condition, now the 
total spin of the three states $*, a = +,0, — , is necessarily St — 1 with S T = +1,0, —1, 
respectively, since all the other antisymmetrized bonds have zero spin. 

These states have again some soliton-like nature, just as in the other approach above, 
e.g., can be explicitly written as 

<T)+ — F Pj-20ij-i d-i pfa-iajJJ /J+l /3 j+1 a j+2 j 3+2 p j+2 a j+3 /io\ 

In each nonzero configuration the subscripts fli and cti+i are antiparallel except for i = j, 
for which f3j = ctj+i =t- Transforming this into the standard S z representation, the nearest 
nonzero spin states on the left and right side of the triplet bond are necessarily + states. 
Otherwise the Nel order of the + and — states is complete on both sides. 

As a variational ansatz, Knabe |H| analyzed a general linear combination of the states, 
Z)jLiCj$", and found that the energy is minimized if Cj = (— For the primary gap of 
the model he obtained an upper bound in the form < A w 0.3571. Correcting a small 
obvious mistake in the numerics of his paper, the correct upper bound is 10/27 ~ 0.3704, 



exactly as in Ref. |L4 . 
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It is quite straightforward to generalize this calculation to arbitrary momentum k. Look- 
ing for translationally invariant trial wave functions, we define 

N 

$ a (Jfe) = Y,j kj ® a j , a=+,0,-. (19) 

3=1 

Knabe's wave function corresponds to k = 7r. With our definition of the (unnormalized) 
trial wave functions a = +,0,—, Knabe's results can be reexpressed in the following 
form 

($•!$•> = I 3 ^(-I)W-i'l+0(l), (20) 

and 

(^\HvBs\n) = 5 J3^ N + 0(1) . (21) 

Using these results, it is now straightforward to calculate the normalization of the transla- 
tionally invariant states $ a (fc) defined by Eq. fll9|) , and the expectation value of the energy 
in these states. In the thermodynamic limit we get 

JV 

($ a (k)\$ a (k)) = e ikU '- j) (^\^,) 

j,3'=0 

1 N 1 9 

= AT V Y e ikr (~) r = / N3 N , (22) 

2 3 5 + 3cos(fc) 

and 

N in 

($ a (k)\H VBS \$ a (k)) = ]T e ik ^{^\H VBS \^) = -N3 N . (23) 



hf=o 

Whence the dispersion of the excitations is again 



(^(k)\H YBS \^(k)) = 25 + 15cos(fc) 
1 ' (<f> a (k)\<S> a (k)} 27 1 J 

This dispersion is exactly the same as that obtained by Arovas et al. in Eq. @. Despite 
the different forms of the wave functions, the identical result for the dispersion relation 
indicates a deep connection between the two approximations. In fact, it is not too difficult 
to show that Sj\Q) can be expressed in a simple form with our configurations: 
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5f|n) = =F(^-^-i), (26) 

therefore differs from $^(A;) in a constant factor only, which is cancelled when the ex- 
pectation value of the energy is taken. 

One may ask the question, which construction of the above two should now be considered 
as the elementary excitation. Eqs. ( p5|) and (|26|) show that Sj\Q) is a simple linear combi- 
nation of the configurations. The revers, however, is not true. cannot be expressed 
with Sj\Q) in a similarly simple way. Therefore, one should conclude that the elementary 
excitations are in fact the moving triplet bonds. 



IV. THE KENNEDY- TASAKI TRANSFORMATION 



A better picture of the above described elementary excitations of the VBS model can be 



obtained by using the nonlocal unitary transformation U of Kennedy and Tasaki [I0|. This 
transforms the antiferromagnetic Hamiltonian into a ferromagnetic-like model and makes 
the Z 2 x Z 2 symmetry breaking in the ground state explicit. First we recall some features of 
this transformation, then show that the above trial wave functions transform under U into 
simple, explicit domain walls. A small inconvenience, arising from the nonlocal character of 
the transformation, is that it can only be used conveniently on chains with open boundary 
conditions. Note, however, that we do not seek for exact solutions but for variational results 
only, so the boundary condition will have no relevance for long enough chains. 



We define the unitary U in the usual way flQ| : Let \s) = \si, «2, • • • , sjv) denote a basis 
state in the S z representation, where = +,0,— stands for the eigenvalues +1,0,-1, 
respectively of Sf. Introducing new variables by 



exp 



(27) 



i=i 

the spin configurations can be given as \s) = |si, S2, • • • , sn)- Note that in \s) all the 0's of 
\s) remain unchanged, while a + or — at site % is flipped or remains unchanged, depending 
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on whether the number of +'s and — 's on sites 1 < I < i is odd or even. The unitary U is 
then defined by 



U\s) = (-l) M ^\s) , 



where M(s) denotes the number of odd sites i on which = 0. 



Let us consider now the VBS Hamiltonian. It can be shown |10| that hj, j 



(28) 



,N-1, 



transforms under U in a relatively simple (local) way into 



hj = UhjU 



-i 
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V i 



-2 



(29) 



1 / 



This, however, does not hold for h^, i-e. for the term that couples the last and the first 
spins of the chain, hjsr = Uh^U~ l cannot be written in a similar form, moreover it does not 
remain local either. This problem can, however, be avoided if we switch to open boundary 
conditions. In this case i^vBS simply reads as 



N-l _ 
3=1 



(30) 



The diagonalization of the above two-site Hamiltonian hj shows that its ground-state 
sector with zero energy is four dimensional, and is spanned by the states </>„ (£></>„, v — 1, 2, 3, 4 
(here no summation is meant over v), where the single-site states are 



(|0)+V2|+))/V3 
(|0> - V2|+))/V3 



(31) 
(32) 
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& = (|0) + V2|-»/V3, (33) 
04 = (|O>-V^|-»/V3. (34) 

Note, that this basis is not orthogonal, since \{4> v >\4>v)\ = 1/3, v' ^ v. 

The ground states of an open chain with iV sites can simply be written as the tensor 
product of the above introduced single-site states, 

* „ = <f>l ® (t>l® ■ ■ ■ <t>u' X ® <Pu , v= 1,2,3, 4. (35) 

These wave functions and the ground-state wave functions described in Eq. (|[) can be easily 
related. On an L-site lattice any linear combinations of the ^'s are ground states. Such 
a linear combination is obtained if U acts directly on the ground state n(a\, of i^vBS- 
For all the configurations in f2(ai,/?z,), the value of a,\ decides unequivocally the sign of the 
first nonzero spin (if a\ =f then it should be a +), and then for a given value of L, (3l fixes 
the parity of the total number of nonzero spins. Using these two observations, the following 
relations can easily be worked out 

UQ(t t) = (-1) 1+L+ ^(3) L/2 ^ 1 - * 2 ) , (36) 
UQ(t |) = (-1) 1+^/21(3)^21^ + ) (37) 

UQ(i, T) = (-l) L+[L/2] (3) L / 2 ^(vI/3 + * 4 ) , (38) 
Un(l, |) = (-l)^ 2 \3) L ^ 3 - * 4 ) , (39) 

where [L/2] denotes the integer part of L/2. 

It can be shown rigorously that the four states in Eq. (|35|) remain the only ground states 
as iV — >• oo, and they converge to four different infinite volume ground states. Note that 
the ground-state degeneracy of the original and the transformed Hamiltonians differ in the 
infinite volume limit. This is a consequence of the nonlocality of U. 

The states in Eq. ( |35| ) have long-range order reflecting the spontaneous breaking of a 
Z2 x Z2 symmetry, the only explicit (local) symmetry of -£/vbs- Introducing the ferromagnetic 
order parameter of the transformed system 
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Og IIO (H)= lim h(S?S«) m , a = x,y,z, (40) 

one obtains that 

= £<<sni = ^ (4i) 

While h( S !}h = +| for *i and *2, we find n( S i)H = ~f for ^3 and *4- Similarly, 
jj(Sf)jj = +| for \l>i and ^3, and a(^i)n = ~ § f° r ^2 and ^4. The appearance of the 
long-range order in the transformed Hamiltonian corresponds to a non-vanishing value of 
the string order parameter in the original system, since by the equivalence 

0: trin JH) = Og IIO (H) for a = x,z, (42) 



the string order transforms into a ferromagnetic order under U ||10|| . This equivalence does 
not hold for the y component. 



V. SOLITONS IN THE KENNEDY-TASAKI TRANSFORMATION 

Our aim now is to show how our soliton configurations a = +, 0, — transform under 
the unitary transformation. It will be found that in the transformed model they are explicit 
domain walls separating regions with different ground states ty u . Since we work now with 
open boundary conditions, the two loose spin-1/2 variables at the left and right chain ends 
should be retained explicitly. 

For this purpose, we cut the chain into £ and 1Z parts as in Section III, and define new 
unitary operators. Since a configuration \s) can be written as an explicit tensor product of 
the left and right states, 

\si, s 2 , ■ ■ ■ , s N ) = \s 1 ,s 2 ,..., Sj) <g> \s j+1 , s j+2 , ■ ■ • , s N ) . (43) 

the action of U on | s) can be given in the form 

U\s u s 2 , • • • , s N ) = Uf\s!, s 2 , . • • , sj) ® Uf\s j+1 , s j+2 , ...,s N ). (44) 
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Here 

Uf\s Xl s 2 ,..., Sj ) = (-1)^%, s 2 , . . . , ^> , (45) 
Uf\s j+1: s j+2 , ...,s N ) = (-l) M ^\s J+1 ,s J+2 , ...,s N ), (46) 

with |s) defined by Eq. (|7]). M c (s) (M^(s)) is the number of odd sites for which s, = 
and % G £ (i E K). Obviously M c (s) + M n {s) = M(s). 

It is seen that in general the effect of [/* on \sj + i, Sj +2 , . . . , sn) cannot be calculated 
without some explicit knowledge of the configuration on C, since e.g. the first nonzero spin 
of 1Z is flipped according to whether the parity of the number of nonzero spins of £ is even 
or odd. Therefore, we define another unitary Vj 1 , that will be independent of £, in the 
following way: 

V^\s j+1 ,s j+2 ,...,s N ) = (-l) K ^\T j+1 ,T j+2 ,...,^ N ) , (47) 

where 

i-l 

Si, i>j + l, (48) 



cxp 



in s i 

l=j+l 



and K-ji(s) is the number of sites i for which i —j is odd and Sj = 0. Note that the definition 
of Vj 1 is nothing else but that of U, if the sites of the chain are relabelled as i — ► i — j. To 
illustrate the above definitions, compare the following examples 

| + + -0) ® Ug | - 00 - + + 0) = -| + + -0} ® | + 00 - - + 0) , 
| + + -0) ® Vf I - 00 - + + 0) = I + + -0) ® | - 00 + + - 0) , 
| + + -0-) ® 17*1 00 - + + 0) = -|+0 + -0-) ® 1 00 - - + 0) , 
| + + -0-) g> Vf|00 - + + 0) = -|+0 + -0-) ® 1 00 - - + 0) . 

It is easy to see that in general Uj 1 can be expressed by Vj 1 as 

Uf = p{s)P(s)V] 1 , (49) 

where p(s) = ±1 is a sign factor and P{s) is either the identity operator or a general spin-flip 
Sj — > — Sj, j + 1 < i < N, on TZ. To be more specific, let us introduce the notation Qc{s) for 
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the number of sites i (i G C) for which s» = ±1, and similarly introduce Qn(s) for % e 1Z. 
For the operator P, we simply get 

{identity if Qr(s) = even , 
(50) 
spin- flip if Qc(s) = odd . 

As for the sign factor p, it is trivially +1 if j — even, while for j = odd it is 

= (_1)AT- J -Q W ( S ) = _(_ 1 )AT-Q W ( S ) ) J = odd; (51) 

since for odd j, K n {s) counts the 0's on even sites and thus M n (s) + Kn(s) is the total 
number of 0's on 1Z. 

Let us consider now the states $j J,Qij+1 (a;i, (3jsr), /3j,aj + i =|, |, where the valence bond 
between sites j and j + 1 is simply removed. From Eqs. (^) and ( f£3j ) 

r<I>; (n,. iy) = UfndauPj) ® Ufn n {a j+1 ,P N ) . (52) 
Then, using Eq. ( pE5| ) we obtain 

U$?' a < +1 (ax,p N ) = Ufn c (a 1 ,/3 j )® P PV^h.(a j+1 ,/3 N ) . (53) 

The expressions for U^Vtc{. a i^j) an d ^-^^^(cKj+i, /3jv) can be read off directly from Eqs. 
(j36H39"l), using the fact that the number of sites in the left part is L = j, while it is L = N — j 



in the right part. Care has to be taken in the proper account of p and P. Remember that 
these depend on the actual configurations \s). However, fixing a.\ and (3j, the parity of Qc(s) 
is uniquely determined for all the possible configurations in flc(aci, /3j). Without any loss 
of generality we will fix the leftmost spin-1/2 variable to a\ =f and suppose that iV=even. 
Then Q£(s)=odd and thus P is a spin-flip [cf. Eq. (|50|)1 if and only if f3j =]. Similarly, the 
parity of Qn{ s ) is uniquely determined by <x, + i and (3^. From Eq. fl5"T|) we easily get 

if j = even , 

P={ (54) 

if j = dd . 



Using the above results and the fact that a spin flip transforms — > ^3, \&2 ^4 and 
vice versa, one straightforwardly obtains the following relations 

f/$j< T (T, T) = (-) J (-3) Ar / 2 (l/4)(M/f - ® (tf* - ) , (55) 

U&fyt T) = -(-3) JV / 2 (l/4)(M/f - *f ) ® (¥* + ¥*) , (56) 

C/^ T (T, T) = (-3)^ 2 (l/4)(*f + *f ) ® (¥* - ¥*) , (57) 

tre^CT, T) = -(-y(-3)^(l/4)(«f + *£) ® (*? + *?) , (58) 

U^\l I) = (-)i(- 3 )*/ 2 (l/4)(*f - ® (¥* + ¥*) , (59) 

C/^(T, |) = -(-3) JV / 2 (l/4)(4f - «f) ® (¥* - ¥*) , (60) 

UQfHl, I) = (-3)^ 2 (l/4)(*f + *f ) ® (¥? + ¥*) , (61) 

l^(T, I) = -(-) J (-3)^ 2 (l/4)(*f + *f ) ® (¥* - ¥*) , (62) 



where the superscript £ (72.) indicates that the wavefunction ^ refers to the left 

(right) part of the chain. 

Let us insert now the triplet bond in place of the missing valence bond, i.e. symmetrize 



with respect to the superscripts of $^" a:,+1 . Recalling Eqs. (P^-fT^D, we find 

tf$t(T, t) = (-) J '(-3) JV/2 (l/4)(*f ® - tff ® - tff ® tf? + *f ® , (63) 

t) = (-3) JV/2 (l/2v / 2)(^ 2 : (8) *f - tff (8) 4?) , (64) 

Ef$j(T, T) = -(-) J (-3) iV/2 (l/4)(4f (g) + (8) + tff ® + *2 ® *4 ) > (65) 

£/$+(T, I) = (-) J (-3) 7V/2 (l/4)(4f ® *J - *f ® *J + *f ® *f - *f ® , (66) 

f/$°(T, |) = (-3) JV / 2 (l/2v / 2)(^ 2 : ® *f + tff (8) , (67) 

f/$7(t, |) = -(-)'(-3)" /2 (l/4)(*f ® + *f ® - tff ® - *f (8) • (68) 



What we obtained is nothing else but a linear combination of the simplest domain walls 
(g> \t^, 1/ 7^ z/, between sites j and j + 1. By virtue of the £77(2) symmetry of 77vbs 
and the open boundary condition, it is possible to consider some linear combinations of the 
above states in order to get the simplest forms on the right hand sides, e.g., 
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f/($+(T, I) + $+(T, T) - $7(t, I) - $7(t, T)) = (-y(-3)^*f ® , (69) 
c/($t(t, i) - $+(t, t) + *j(t, I) - *J(t, T)) = (-y(-3)^*f ® , (70) 

r/($;(t, i) - *j(t, t)) = (l/v^x-s^f g, ^ . (71) 

Three other similar linear combinations can be composed with ^ (g) \E^, z/ = 1,3, 4, on 
the right side. These final forms clearly demonstrates the solitonic nature of our trial wave 
functions. Note that in accordance with the three degrees of freedom of such a spin-1 soliton, 
there are three kinds of domain walls. For \l/f ® (g) \I/^) only fj{S?)jj changes sign 

as we move from the left region to the right, jj(Sf) # does not change. The situation is just 
the opposite for *f <g> tyf ® Wf). Here only g(Sf) g flips. Then for \Pf <g> *J <g> Hff) 
the expectation values of the magnetization in both directions change sign. 

Now it is easy to see how <3>° [or SjQ, recalling Eq. (|25|) 1 destroys the hidden order Og tring 
(and by symmetry Og tring ) which was anticipated in Sec. III. For definiteness, we fix the 
boundary spins a.\ =|, (3^ =T (other choices can be worked out similarly) and consider the 
expectation value 

(*;(t,t)KJ$;(u)) (79] 

In the case when j < n or j > m, i.e. n and m are in the same domain, the domain wall 
has no effect and the expectation value is 4/9. On the other hand when n < j < m, it is 
straightforward to obtain, using the Kennedy- Tasaki transformation, the following result in 
the thermodynamic limit 

($°(U)KJ$°(U)) _ (^°(T,T)|^|^°(T,T)) 
ll^(T,T)ll 2 ll^?(T,T)ll 2 





5 *f - *f € 




- *f <g 


4 






(g) *f - *f g) ^ll 2 
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where we used Eq. (|64] ) and the asymptotic orthogonality of the different ground states 
| (^f u '\^f u ) | — > (v' 7^ z/) if iV —>■ oo. In fact, the presence of the domain wall flips the 
expectation value of cr^ m . 

Finally we show that in this formalism the dispersion relation of Eq. ( |T0"D can be obtained 
in a very elegant way. We can start from e.g. the trial wave function 
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N-1 

\k) = £ e ikj \j) > (74) 

3=1 

with \j) = (f)\ ® ® . . . (f)[ ® (f) J 2 +1 <E> 4>2 +2 ® • • ■ 02^- m this form /c is a variational parameter 
rather then a momentum, because of the open boundary condition. In the thermodynamic 
limit, however, the boundary condition should not matter (although it might bring a constant 
momentum shift q in the final result, since we have the freedom to redefine |j) with e.g. an 
arbitrary phase factor \ j) — > e iq: >\j)), and the variational energy as a function of k is in fact 
the dispersion of the excitations in this simplest domain wall approach. 

The computation would proceed similarly to that in the previous section. As we have 
seen there, the important quantities are (j\f) and (j\HvBs\j')- However, in this case they 
are trivial because of the tensor product form. A straightforward calculation gives 

W) = (-l) lj ~ fl , (75) 

where we used that (fa\fa) = (fa\fa) = 1 and (4>i\<p2) = ~ 1/3, and in a similarly simple 
way 

20 

0'1-^vbsIj") = <W(0i ® fa\hj\4>i ® fa) = > (76) 

which is easily obtained from the explicit form of the two-site Hamiltonian hj. Apart from 
a factor of 3 N /2 these results are identical to those in Eqs. (|20|) and (|21|). Therefore, they 
also lead to the same dispersion e(k). 

VI. CONCLUSION 

In summary, we studied the elementary excitations in the valence-bond point of the 
5=1 bilinear-biquadratic model. Numerical calculations on finite-size systems were used 
to predict the spectrum in the thermodynamic limit. The lowest-lying excited states above 
the k = singlet ground state form a discrete triplet branch with a minimum at k = tt. 
Near this minimum this branch is separated from the higher-lying scattering continuum. 
The energy needed to excite the lowest k = excitation was found to be twice the gap value 
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at k = it. Similarly, the energy of the next lowest excitation at k = it is three times the gap 
value. These excitations belong to the continuum and can be interpreted as being composed 
of two or three S = 1 elementary excitations. 

Comparison with the numerical results show that the separate branch of excitations can 
be reasonably described with a trial wave function, where one singlet bond is replaced by 
a moving triplet bond. In the representation where the configurations are given in terms 
of the S z eigenstates of the spins, a triplet bond in the sea of singlet bonds has a solitonic 
character. In the dilute system of + and — spin states there is a single domain wall. While 
this feature is hidden in the usual valence-bond description, it becomes apparent when the 
nonlocal Kennedy- Tasaki transformation is used. We have shown that the approximate 
wave functions of the excited states transform into explicit domain walls in the transformed 
system. 

This research was supported in part by the Hungarian Research Fund (OTKA) Grant 
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FIGURES 

FIG. 1. Low-lying eigenvalues of the VBS model plotted vs momentum k, for a chain with 
N = 16 sites. Labels denote the total spin St of the states. Dashed line shows the energy of the 
trial state with the moving hidden soliton. 

FIG. 2. Energy gaps A^, As/2, and Ac/3 plotted vs 1/N. Dashed lines indicate the suggested 
large N behavior. □ shows the energy of the trial wave function at k = it. 
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